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The formalism of Ursell operators provides a self- consistent integral equation 
for the one-particle reduced operator. In three dimensions this technique 
yields values of the shift in the Bose-Einstein condensation (BEC) transition 
temperature, as a function of the scattering length, that are in good agreement 
with those of Green's function and quantum Monte Carlo methods. We have 
applied the same equations to a uniform two-dimensional system and find 
that, as we alter the chemical potential, an instability develops so that the 
self- consistent equations no longer have a solution. This instability, which 
seems to indicate that interactions restore a transition, occurs at a non-zero 
value of an effective chemical potential. The non-linear equations are limited 
to temperatures greater than or equal to Tc, so that they do not indicate 
the nature of the new stable state, but we speculate concerning whether it 
is a Kosterlitz-Thouless state or a "smeared" BEC, which might avoid any 
violation of the Hohenberg theorem, as described in an accompanying paper. 
PACS numbers: 03.75. Fi,05.30.Jp,05.70.Fh,67.40.Db, 68. 35.Rh. 



1. INTRODUCTION 

Interest in Bose systems and Bose-Einstein condensation (BEC) has 
been expanding in recent years, stimulated in part by the experimental work 
on alkali atoms in magnetic traps .□ Much of the theoretical work has involved 
mean-field theories, such as those based on the Gross-Pitaevskii equationB 
However, some theoretical analyses have recently reconsidered the proper- 
ties of three-dimensional (3D) homogeneous systems with the inclusion of 
dynamic correlations that go beyond mean-field theory. Here we extend the 
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methods to dilute systems in two-dimensions (2D). 

In 3D, these techniques have provided computations of the change in the 
BEC transition temperature as a function of the repulsive interactions in a 
dilute Bose gas.B'a'Q Ref. a reports results from an Ursell operator approach,Q 
which are confirmed by non-perturbative Green's function methods-u'El The 
two techniques give equivalent non-linear self-consistent equations for the 
self-energies and the distribution functions. Each method finds that a dilute 
Bose system with repulsive interactions reduces its self-energy by spatial 
rearrangements of particles with low velocities, and that the \ow-k quasi- 
particle energies behave as k'^~^ near the transition. (The constant t] is 
equal to 1/2 in lowest order approximation, and smaller in higher order.) 
Both methods predict that there is an increase of the critical temperature 
as the repulsive interactions are turned on, a result confirmed by quantum 
Monte Carlo calculations .13 

There may soon be a number of experiments involving 2D systems for, 
say, alkali atoms or for spin-polarized hydrogen in traps or in homogeneous 
settings. The usefulness of the Ursell approach in 3D suggests studying 2D 
systems by the same approach. Here we solve the self-consistent equation 
for the self-energy as a function of the chemical potential. In 3D, as the 
chemical potential is increased from large negative values, one reaches a 
critical value where the phase transition takes place; an effective chemical 
potential that includes the effects of mean-field and dynamic correlations 
vanishes signaling the transition. We carry out a similar program here and 
find that at a finite value of the effective chemical potential the equations no 
longer have a self-consistent solution, again signaling some sort of transition. 

The Hohenberg theorems states that a 2D system cannot undergo a 
BEC into a single state, but recent work shows that a "smeared" conden- 
sation is not forbidden by this theorem.113 In this mode of BEC, no single 
state is extensively occupied, but each of the states in a band has a large, 
although non-extensive, occupation, with enough states in the band that 
its total occupation is extensive. Furthermore, oneknows that superfluidity 
occurs in 2D via a Kosterlitz-Thouless transition.Ei iLis even possible that 
the KT transition and a smeared BEC are equivalent .O Since our theory has 
not yet been extended to regimes in which there is macroscopic occupation 
of a single state or of a band of states, we are unable to tell what transition 
is signaled by our instability. 
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2. URSELL OPERATOR APPROACH 

We use the method of Urseh operators,^ an approach that allows the 
generalization of the usual cluster expansion to all orders in statistical effects, 
while still allowing the possiblity of including interactions in lower order. A 
basic quantity is the second Ursell operator given by 

[/2(1, 2) = e-^^^(i'2) - e-/5{^i(i)+^i{2)) (1) 

where Hi is a single-particle Hamiltonian and H2 a two-particle Hamiltonian. 
A particular problem in 2D is that the matrix elements of this operator do not 
approach a constant at small k; indeed, for hard-core scattering, the leading 
term is proportional to l/ln(A;a), where a is the hard-core radius. However, 
we will see that we can, to good approximation, avoid that complication and 
use a k-independent U2 matrix element, as in 3D. 

The Ursell operator method gives an integral equation for the one- 
particle reduced density pi with Fourier Ijiansform pk- When the occupation 
of all states is small (no BEC), one finds,0 to lowest order in the interaction, 

that pk is given by a free-particle distribution function with a self-consistent, 

-1 
, where 

The lead- 



Pisk-fi-k) _ \ 



fe-dependent effective chemical potential: pk 

£k = h k"^ /2m and flk = l^—^P-k is the effective chemical potential 

ing correction to the chemical potential is given by PApk = — In [1 — A.^^] , 

where 

^^^ = -(2^y^ ^n(ik-k1a) ^') 

with a equal to a times a constant of order unity. (In 3D the log factor does 
not occur and so A^ is independent of k.) We have solved these equations by 
iteration and find that A^k is a very weak function of fc, so that one might as 
well set the logarithm to a constant. A constant correction to the chemical 
potential or, equivalently, the single-particle energy is simply a mean-field, 
which can cause no fundamental change, and certainly no phase transition, 
in the Bose gas. Ignoring the logarithm is equivalent to using a momemtum- 
independent matrix element of U2, which is what we use in the next order 
approximation. 

If we add corrections second order in the interaction we necessarily 
develop velocity-dependent terms in the distribution function. We find that 
the effective chemical potential change is 

f3A^ik = -ln[l-ACo-^Ck] (3) 

where, with our assumption of constant interaction, we now have 



ACo = j^ I dkpke^^^'= (4) 
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with a a parameter measuring the strength of the interaction. The second- 
order correction in Eq.(^) is 

The quantities A£q and (5^k a-re equivalent to self-energies in the Green's 
function approach .□ 

There are a variety of calculational procedures that can be followed to 
compute (^^k- If one assumes that p^ is dependent only on the magnitude 
of k, then this quantity can be reduced to a four-fold numerical integral. 
By changes of variables this integral can be done in several ways. Alterna- 
tively, one can introduce the Fourier transform of PkC^ ; since Eq.(g) is a 
convolution, the computation reduces to doing the Fourier transform and a 
one-dimensional integral. The procedure used is to pick a value of the real 
chemical potential, compute the interaction corrections, and iterate. As in 
the 3D case, 5^k has a minimum at fc = and rises monotonically. As one 
increases /x from very negative values, as would happen in 3D in a system 
approaching criticality, the minimum in (5^k deepens and the curvature at 
A; = sharpens. The value of \flk\ decreases during this process. This feature 
is caused by dynamic correlations that lead to rearrangements of particles 
with small k. 

Unexpectedly, at a particular value of /i, with \jlk\ still nonvanishing, 
the iteration becomes unstable and the self-consistent solution disappears 
as shown in Fig. 1. This result is quite robust and occurs independently of 
the numerical method used. At the instability \6^o\, instead of converging 
to a self-consistent value as one iterates, increases rapidly until the process 
predicts a positive value of pQ, at which point the iterations must halt. 
Passing through //q = would seem to indicate the onset of BEG. Since the 
theory has not been extended into any condensed regime we are unable to 
say what the new stable phase might be. 

Since one knows that there is a KT phase in 2D, one might expect 
that our theory is showing an instability towards that phase. A second 
possibility arises from the proof in an accompanying paper ,t3 which shows 
that the Hohenberg theorem does not rule out what we term as a "smeared" 
BEG, as defined in the Introduction. It will be interesting to extend the 
theory to cover such a possibility. 

One can go to higher orders in the approximation scheme and sum all 
"bubble" diagrams, in which case one finds 

'^^ - 2(2vr)4 y ^^ ^-^ ^ j 1 _ ^ J dqp,Pk+.+lc'e^(^'^<'+^-^...-) J 

(6) 
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We find that the instabihty persists, although the inter ations now have a 
somewhat different behavior. As one iterates, values of A^o and J^k conspire 
to make the subsequent values of these quantities oscillate about an average 
value. Either this oscillation dies out eventually to a self-consistent value or 
grows until it becomes completely unstable. 



3. DISCUSSION 

Because of the success of the Ursell approach in 3D, namely, its good 
agreement with the results of other powerful approaches, we have confidence 
that it should yield reliable results in 2D. Our results show an instability in 
this case, but we are unable, at this time, to determine if this is the signal 
of a true phase transition, what is the order of such a transition, or what 
might be the stable state to which the system proceeds. It is known that 
perturbation theory breaks down near the BEC phase transitiona lI, and 
that applies here as well. Nevertheless, the unknown sum of the diagrams 
associated with interparticle correlations must give some negative result, 
because the mean-field theory can be seen as a variational evaluation with 
the same wave function as an ideal gas, and the inclusion of correlations can 
only reduce the free energy and bring the system closer to intability. As /x 
and the density increase, it seems then that an instability must ultimately 
take place. Possiblities are that we are seeing a KT transition or, perhaps 
more interestingly, a transition to a "smeared" Bose-Einstein condensation, 
which our accompanying article shows is possible.t3 Thus the present paper 
is really a progress report and we will need further research to resolve the 
questions it raises. 
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5. FIGURE CAPTION 

FIG. 1 Change in the effective chemical potential for fc = as a func- 
tion of the real chemical potential. For the interaction strength used the 
instability occurs at f3^ > 0.13, beyond which the plot shows successive final 
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iteration values just before failure rather than self-consistent values. BEC 
criticality occurs when fi + A^ = 0, i.e., when the A/Uq touches the diagonal 
H line. 
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